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The deformation of two-dimensional vortex patches in the vicinity of fluid boundaries is investi- 
gated. The presence of a boundary causes an initially circular patch of uniform vorticity to deform. 
Sufficiently far away from the boundary, the deformed shape is well approximated by an ellipse. 
This leading order elliptical deformation is investigated via the elliptic moment model of Melander, 
Zabusky & Styczek [M. V. Melander, N. J. Zabusky & A. S. Styczek, J. Fluid. Mech., 167, 95 
(1986)]. When the boundary is straight, the centre of the elliptic patch remains at a constant 
distance from the boundary, and the motion is integrable. Furthermore, since the straining flow 
acting on the patch is constant in time, the problem is that of an elliptic vortex patch in constant 
strain, which was analysed by Kida [S. Kida, J. Phys. Soc. Japan, 50, 3517 (1981)]. For more 
complicated boundary shapes, such as a square corner, the motion is no longer integrable. Instead, 
there is an adiabatic invariant for the motion. This adiabatic invariant arises due to the separation 
in times scales between the relatively rapid time scale associated with the rotation of the patch 
and the slower time scale associated with the self-advection of the patch along the boundary. The 
interaction of a vortex patch with a circular island is also considered. Without a background flow, 
conservation of angular impulse implies that the motion is again integrable. The addition of an 
irrotational flow past the island can drive the patch towards the boundary, leading to the possibility 
of large deformations and breakup. 



I. INTRODUCTION 

Many problems in geophysical fluid dynamics feature quasi-two-dimensional flows due to the presence of 
rotation, stratiflcation, or shallow water; and a ubiquitous feature of such flows is the presence of large-scale 
vorticity structures. Two examples can be found in the Atlantic Ocean: an instability of the Gulf Stream 
leads to the formation of Gulf Stream rings of relatively hot or cold water—; and the flow of water out of 
the Mediterranean Sea leads to the formation of Meddles, which are rotating regions of relatively warm 
and salty wateriSi^. In both cases, the vorticity structures can have life times of more than a year, and are 
important in the transport of heat and salt across the Atlantic. During their life times, these structures can 
interact with both coastal and submarine topography, which leads to deformation of the structures. For 
example, Richard & Tychensky^ observed the interaction of Meddies with underwater seamounts, and found 
that such interactions could cause the Meddy to break up. 

The simplest formulation within which to study such flows is the two-dimensional incompressible Euler 
equations. This formulation has been used by Johnson & McDonald^i^ who modelled vortex patches as point 
vortices in order to investigate the paths of ocean vortex patches in the vicinity of boundaries such as islands 
or coastal gaps. However, with the exception of a few contour dynamics calculations (e.g. Appendix C of 
Thompson^i and Crowdy & Surana-), there has been little focus on the deformation of vorticity structures 
caused by the presence of boundaries. This motivates the more systematic investigation presented here. 

As a model for more complex vorticity structures, we consider (initially circular) patches of uniform 
vorticity and utilise the elliptic model of Melander, Zabusky & Styczek-. This model has been widely used, 
with successful applications including the study of symmetric vortex merging both with^ and withoulji^ an 
imposed background shear flow. We favour the simplicity of this low-order model over the potential for 
increased realism from models with varying vorticity such as the elliptical model of Legras & Dritschel'^i'. 
When more complex vorticity structures have been considered, the action of an external straining flow 
has been found to strip off outer layers of vorticity leading to the steepening of vorticity gradients and 
consequently to structures that more closely resemble the uniform patches of vorticity that we consider 
herc^^. 

In Section |TT1 we describe the elliptic model of Melander, Zabusky & Styczek^, and we also give a brief 
description of the contour dynamics method that can be used to calculate numerically the full evolution of 
patches of uniform vorticity. Then, in Section IIIIl we apply the elliptic model to the motion of an initially 
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circular vortex patch in the vicinity of a straight boundary. The resulting motion is integrable, and is shown 
to reduce to the problem of an elliptic vortex patch in constant strain, which has been analysed by KidaF^. 
Comparison with contour dynamics solutions confirms that the elliptic model is a good approximation to the 
motion, even when the separation of the vortex patch from the boundary is comparable to the dimensions 
of the vortex patch. In Section llVl we consider the motion of a vortex patch around a corner, which is no 
longer integrable. The time scale associated with the rotation of the patch is found to be short compared 
to the time scale on which the straining flow experienced by the patch varies. This separation of time 
scales in a Hamiltonian system implies the existence of an adiabatic invariant, which is an (approximate) 
constant of the motionii. We briefly consider the effect of adding a background straining flow in Section |Vl 
This flow can trap vortex patches in the corner. We use Poincare sections to reveal the underlying chaotic 
nature of the elliptic model. Then, in Section IVH we consider the deformation of a vortex patch due to a 
circular island. Without a background flow, the motion around an island is found to be integrable. When 
a background flow is added, significant deformation, and possible breakup of vortex patches is observed as 
they are driven towards the boundary. 



II. METHODS 
A. Elliptic model 

A vortex patch of uniform vorticity, in the vicinity of a boundary, will generate a flow field u(x) which 
can be decomposed as u(x) = Uoo(x) + ub(x), where Uoo(x) is the flow fleld that the patch would generate 
in an unbounded geometry, and ub(x) is the perturbation to that flow due to the presence of the boundary. 
For an elliptical vortex patch, Uoo(x) simply causes the patch to rotate while remaining the same shape^^. 
The flow due to the boundary, ub (x) , can be decomposed about the centre of the vortex patch as a uniform 
advection, a straining flow, and higher order corrections. The uniform advection simply translates the vortex 
patch with no change of shape. The straining flow will change the shape of the vortex patch, but an intially 
elliptical vortex patch will always remain elliptical in shape^^. Thus, provided the higher order corrections 
can be safely neglected, the effect of a boundary on an initially circular vortex patch is to deform it into 
an elliptical shape. This observation motivates the use of an elliptic moment model to approximate the 
evolution of such a vortex patch. 

When the boundary shape is sufficiently simple, the motion can be solved via the method of images. 
In this case, the motion is simply that of a collection of elliptical vortex patches. A Hamiltonian elliptic 
moment model for such a collection of vortex patches was developed by Melander, Zabusky & Styczek*. In 
Section [II CI we will recap the equations of motion associated with this model, but first, we briefly consider 
how to represent an elliptic shape mathematically. 



B. Representing an ellipse 

For an ellipse of area A, two further independent quantities are needed to describe the shape. Perhaps 
the simplest and most physically intuitive pairing is {r,9), where r > 1 is the ratio of the length of the 
major axis of the ellipse to that of its minor axis, and 9 is the angle that the major axis makes to some 
predetermined direction (here we will always measure angles relative to the positive a;- axis). 

A second representation can be obtained by considering the quadratic moment of an elliptic patch 



where Xc is the centroid of the ellipse and the integral is taken over the area of the ellipse. This can of 
course be written in terms of r and 9 



Since there are three quantities qxxi Qxy and qyy representing the shape, they can not all be independent; 
they satisfy the relation det(q) = A'^/{Att). 




(1) 




(2) 
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To facilitate a compact representation of the elliptic model in the next section, we also define a matrix p, 
closely related to q, as 

f q^y {<lyy-<lxx)l'2\_^, i.( sm{2e) - cos(26')\ , . 

P - \{,qyy - gxx)/2 -q^y ) Svr ^ ^ V- cos(20) - sin(2^) ) ' ^""^ 



C. Model equations 

The motion of well-separated elliptical vortex patches is approximated by the Hamiltonian elliptic model 
of Melander, Zabusky & Styczek^. This model approximates the evolution of a collection of N vortex patches 
with vorticity Wi, area Ai, centroid and elliptic shape (r^, 9i). The associated Hamiltonian is given by the 
excess energy of the systemic 

N 

where ipip^) is a stream function for the flow. Each integral over a vortex patch can be evaluated as 

H. - ^ 1" (^) + 1 '"(««) - 1 W 



self 

where 



point elliptic 



_ ujjA, ( sin(2%) -cos(26ly)\ , , 

" 2^Ef^A-cos(2%) -sin(2%)y W 

is the rate of strain at the centre of patch i due to the leading-order point vortex contribution from patch 
j, Rij = |xi — Xjl is the distance between the centres of patches i and j, and 9ij is the direction of x^ — Xj. 
Terms of order and higher are neglected in the calculation of the Hamiltonian. This approximation is 
valid provided Rij is large compared to the longest dimension of patches i and j. 

The Hamiltonian consists of three parts: a 'self energy, a 'point' energy and an 'elliptic' energy. The self 
energy is associated with the self-induced rotation of the vortex patches. The point energy is associated 
with the leading-order motion of the patch centroids due to point vortex interactions. The elliptic energy is 
associated with changes in shape due to strain from other patches and also with the next-order corrections 
to the motion of the patch centroids. We choose to write the 'elliptic' contribution to the Hamiltonian in 
a slightly different form to that given by Melander, Zabusky & Styczek. Our form has the advantage of 
emphasising the role of strain in changing the shape of the ellipse. 

The model is completed by the evolution equations 



1 dH 
uiiAi dxi 
1 dH 

uJiA^ dui 
8tt r? dH 
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,A^ 1 - r? 89, " 



(7a) 
(7b) 
(7c) 
(7d) 



Melander, Zabusky & Styczek derive these equations by inspection from their approximation of the Eu- 
ler equations. Alternatively, they can be derived directly from the Hamiltonian formulation of the Euler 
equations (see Meacham, Morrison & Flierli^). 
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We note that the model can be easily extended to include the effect of an irrotational background flow. 
The additional excess energy due to an irrotational background flow with stream function is given by 

N 

Hb^-Y,^^ / ^(x)d'x. (8) 

^=l •^'4. 

The elliptic approximation is then obtained by evaluating these integrals and neglecting any moments that 
are higher than second order. 

A vital assumption in the use of this Hamiltonian model is that the effects of viscosity can be neglected. 
Viscosity causes the patch vorticity to diffuse outwards and leads to the formation of viscous boundary 
layers along fluid boundaries; but perhaps most critically, it breaks the Hamiltonian nature of the system. 
Such effects can be neglected, and the system treated as Hamiltonian, provided the appropriate Reynolds 
numbers are large. A patch of vorticity uj and typical dimension R generates a typical velocity ujR. Thus 
diffusion of the patch viscosity can be neglected provided the Reynolds number ujR^ /v is large. The effects 
of viscosity at fluid boundaries can also be neglected provided another Reynolds number luRS/ ly, based on 
the separation S between the patch and the boundary, is also large. 



D. Contour dynamics 

The elliptic model is only valid provided the vortex patches, or equivalently, the patch and boundary, are 
sufficiently far apart. In reality, this separation does not have to be very large for the model to be a very 
good approximation. Nevertheless, it is informative to compare the model against the full solution. This full 
solution can be calculated relatively efficiently via the method of contour dynamics^'* as briefly described 
below. 

The stream function due to a patch of uniform vorticity is given by an integral, over the area of the patch, 
of the point vortex Green's function 

V'(x) = -|^^ln(|x-i|)d2i. (9) 

To calculate the fluid velocity due to this patch, derivatives of the stream function must be evaluated. In 
the contour dynamics method, the divergence theorem is used to convert the expression for these derivatives 
from an integral over the area of the patch to an integral around the boundary of the patch: 

VV(x) = ^^Viln(|x-x|)d2x 

f ln(|x-x|)dn. (10) 

JdA 

Thus the velocity of points on the vortex patch boundaries, which are required to evolve the patch positions, 
can be calculated by integrating around only the boundaries of the vortex patches. 



III. STRAIGHT BOUNDARY 



The simplest interaction of a vortex patch with a boundary is the interaction with a straight boundary. 
The motion of a vortex patch in the half-plane y > above a boundary at ?/ = is equivalent, via the 
method of images, to the motion of a vortex patch lying in ?/ > combined with its mirror image in y = 0. 
This image lies in the half-plane y < and has the opposite sense of vorticity (Figure [Ij . The Hamiltonian 
for this system can be constructed using Equation ([5]) . The resulting Hamiltonian is eight-dimensional with 
four parameters describing each of the original vortex patch and its image. The Hamiltonian governing only 
the motion of the original patch can be obtained by writing the parameters of the image patch in terms of 
those of the original patch, and by considering only the excess energy associated with the half-plane y > 0, 
which is simply half the excess energy of the full plane. Thus the Hamiltonian for the original patch is 

Hiy.r^O) . M)!ln ((i±^) - (^ln(2,.) - |p(.,.): E(,.), (11) 
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FIG. 1. The motion of a vortex patch above a straight boundary is equivalent, via the method of images, to that of 
a vortex patch and an image vortex patch, of opposite vorticity, obtained by a reflection in the boundary. 

where 

Due to the translational symmetry of the geometry, the Hamiltonian is independent of Xc- This symmetry 
imphes the existence of an additional conserved quantity, which is the distance yc from the boundary of 
the centroid of the ellipse. For the four dimensional elliptic model, the presence of an additional conserved 
quantity implies that the motion is integrable. Furthermore, since yc is constant, the evolution of the elliptic 
shape is simply that of an elliptic vortex patch in a constant straining background flow. The steady solutions 
for an elliptic vortex patch in constant strain were found by Moore & Saffman^^. The unsteady case, which 
is of interest here, was later analysed by KidaJ^ (who actually analysed the more general problem of an 
ellipse in shear flow). We can recover his result for the subcase where there is no background vorticity by 
substituting in for p and E in Equation ([TT|) to obtain 

e r2 - 1 \ Asr 

where e = {Acu) / (Siry'^) is the strength of the straining flow, and s is some function of the constants H and 

Vc- 

Kida found that there are three fundamental types of motion possible: indefinite extension, in which 
the long axis of the ellipse aligns itself with the axis of extension of the strain {0 = — 7r/4) and r ^ oo; 
rotation, in which r is bounded and 6 varies monotonically; and nutation, in which r is bounded and 6 
oscillates about 9^0. For sufficiently strong strain, e/w > 0.150, the only possibility is indefinite extension. 
For 0.123 < e/oj < 0.150 both extension and nutation are possible. Then for sufficiently weak strain, 
e/oj < 0.123, all three types of motion are possible. Figure [2] shows possible motions in the r-9 plane when 
e/uj = 1/18 — 0.056. The closed contour represents a nutating motion. The contours which arc monotonic 
in 6 represent rotating motions. In between these two types of motion lies the motion of an initially circular 
patch. These three possibilities were shown in Meacham et al,'^^ to be conveniently determined by the 
intersection of a Casimir hyperboloid with the energy surface, just as the motion of the free rigid body is 
determined by the intersection of the angular momentum sphere with the energy ellipsoid. 

A quantitative measure of how much a boundary deforms vortex patches is the maximum deformation, as 
measured by r, experienced by an intially circular vortex patch. An initially circular vortex patch, of radius 
R, will be indefinitely extended if yc/R < 1-01 (corresponding to e/ut > 0.123) or return to its original 
circular shape in a periodic motion if yc/R > 1.01. For the latter case. Figure [3] shows the maximum 
deformation rmax of the intially circular vortex patch in terms of the distance of its centroid from the 
boundary. The maximum deformation falls off rapidly with distance from the boundary. This fall-off is a 
consequence of the 1 /y^ decay of the straining flow acting on the patch as yc increases. 

A. Model validation 

The elliptic model is only valid if the vortex patch and its image are well separated, but in reality, it 
provides a good approximation to the motion of an initially circular vortex patch even when it is remarkably 
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FIG. 2. Possible motions in the r-9 plane as predicted by the elliptic model with e/ui = 1/18, corresponding to 
Uc/R ~ 1.5. The contour that intersects r = 1, at which point 9 is undefined, represents the motion of an initially 
circular patch (solid line). Any contour within this circular motion contour represents a nutating motion (dot-dashed 
line). Contours outside the circular motion contour represent vortex patches undergoing rotation (dashed lines). For 
sufHciently large values of r, larger than those shown here, vortex patches are extended indefinitely. 



close to the boundary. Figure H] shows the evolution of such a vortex patch with j/c = 1.5i? both under the 
eUiptic model and the full contour dynamics solution. There is very good agreement between the two, with 
the most notable difference being a slight phase drift with time. 

Even for yc — R, when the contour dynamics solution shows that the vortex patch is stretched out into a 
non-elliptical shape, the elliptic model gives a remarkably good approximation to the overall region occupied 
by the vortex patch at short times (Figure O . 



IV. MOTION AROUND A CORNER 

Since the motion above a straight boundary is periodic, the only way that large deformations to the 
vortex patch can occur is if the separation of the patch from the boundary is comparable to its radius. For 
more complicated boundaries, we might expect deformations to the patch shape to accumulate over time. 
In which case, even vortex patches whose separation from the boundary is large compared to their radius 
could eventually accumulate large deformations. This motivates us to consider the motion of a vortex patch 
in a quarter plane (a; > and y > 0). 

Since a quarter plane does not possess the translational symmetry of the half plane, the motion of a vortex 
patch within a quarter plane is unlikely to be integrable under the elliptic model. This could be determined 
by a Melnikov analysis, although such analyses in the present context can have special difficulties. (See 
Ngan, Meacham & Morrison^ and references therein.) Nevertheless, the geometry is still simple enough to 
allow the use of the method of images. In this case, two reflections must be made: one in y = 0, and one in 
X = (Figure [H]). This gives three image vortex patches, two of which have opposite vorticity to the original 
patch. The Hamiltonian for the motion of the original patch (a quarter of the excess energy for the full 
plane with images) is given by 

(-£±1)!) _ £±11! i„ ( .py^ _ " (, „j, (14) 
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FIG. 3. Maximum deformation rmax of an initially circular patch of vorticity in terms of the distance of its centroid 
from the boundary, as given by the elliptic model. The model predicts that the patch is extended indefinitely when 
y,/R < 1.01. 
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FIG. 4. Motion of an initially circular patch with centroid at yc = 1.5R. Patch motion is from left to right at time 
intervals of 5/uj. Contour dynamics solution is given by solid lines, and the elliptic approximation by dashed lines. 



where the strain rate 



E(a;,y) 




y- 



2a: 1/ 



(15) 



now depends on both Xc and yc- 

A patch of positive vorticity will propagate from the region y ^ x, around the corner, and then towards 
the region x ^ y; and vice- versa for a patch of negative vorticity. When either Xc ^ j/c or y^ ^ Xc the 
evolution of the patch is well approximated by that of a patch above a straight boundary, for which the 
distance of the centroid from the boundary is fixed. We define yoo and Xoo to be the fimits of j/c and Xc as 

— >■ oo and yc ^ oo respectively. 

Figure [7] shows the evolution of an initially circular vortex patch with Xoo — 1.5i? moving around the 
corner. The most notable feature of this evolution is the similarity of the periodic variation of r with t 
before and after the interaction with the corner region. The motion afterwards is almost exactly the same, 
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FIG. 5. Even when the elhptic model is not vahd, it still gives a good approximation to the region occupied by the 
vortex patch at short times. Here we show the contour dynamics and elliptic model solutions with yc — R at time 
intervals of 
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FIG. 6. The motion of a vortex patch around a corner is equivalent, via the method of images, to that of a vortex 
patch and its image vortex patches obtained by reflection in the boundaries j: = and y — 0. The two image vortex 
patches obtained from a single reflection have the opposite vorticity to the original patch, the third image patch has 
the same vorticity. 



periodically circular, motion that the vortex patch was executing before it approached the corner. This 
similarity suggests that Xoo ~ Hoc, and indeed we find that \xoo — 2/oo| = 4 x 10~^R. The approximate 
equality between Xoo and i/oo, and the corresponding similar periodic evolution of r, are generic. It appears 
to occur for all initially circular patches with x,^ > l.Oli?, i.e. all initially circular patches that are not 
indefinitely extended when at a similar separation from a straight boundary (Figure [3]) . The explanation 
for this phenomenon lies in the existence of an adiabatic invariant. 




FIG. 7. Evolution of an intially circular vortex patch in a quarter plane with initial centroid (1.57?, lOOiZ). The time 
origin f = is defined to occur when Xc = yc- (a) Path of the patch centroid. (b) Variation of elliptic shape with 
time. Far away from the corner, the patch undergoes a periodic motion associated with the evolution of a patch 
above a straight boundary. The periodic evolution of the patch shape before and after the corner is almost exactly 
the same. 

A. Adiabatic invariance 

Adiabatic invariance occurs for two-dimensional Hamiltonian systems within a slowly varying environment. 
Due to the separation of time scales between the motion of the system and the variation of its environment 
there exists an (approximate) invariant of the motion referred to as the adiabatic invariant^*. Strictly 
speaking, the adiabatic invariant is not an invariant but rather a quantity which varies by only a very small 
amount. It is constant to all orders in e, where e is defined to be the small ratio of the time scale of the 
system to the time scale on which the environment varies, but can vary by an exponentially small amount. 

Returning to an elliptic vortex patch in a quarter plane, we observe that there are two time scales 
associated with the motion of the vortex around the corner: the time scale on which the patch rotates, Tj-ot! 
and the time time scale on which the strain experienced by the patch varies, Tstrain- From Figure [ZUb), 
we see that Trot is the shorter of these two time scales. The variation in strain is largely confined to the 
period —40 < cot < AO but during this period the patch rotates multiple times. This separation of time 
scales, while not an order of magnitude, is still sufficient to lead to the phenomenon of adiabatic invariance. 
(That an elliptic vortex patch in a slowly varying strain can behave adiabatically was first observed by 
Legras & Dritschel^^.) 

The separation of time scales allows us to consider the evolution of the vortex patch shape as a two- 
dimensional dynamical system with a slowly varying parameter instead of as part of the full four-dimensional 
elliptic model. This slowly varying parameter is the strain rate E experienced by the vortex patch. The 
Hamiltonian for the two-dimensional system is given by 

i/,h,pe(r, e, t) ^ In - |p(r, 9): E(ec.i). (16) 

We write E — E(ewi), where e ^ 1, to emphasise that the variation in strain rate is slow compared to the 
rotation of the patch. The theory of adiabatic invariance applies to this system and implies the existence of 
an adiabatic invariant that is constant to all orders in e. Furthermore, when E is constant, the value of this 
adiabatic invariant is simply given by the action 



/(P,'?,E) = 




where 
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FIG. 8. (a) Evolution of the action / for the motion shown in Figure [T] The values of I before and after the corner 
interaction are almost exactly the same, (b) The difference in action A7 before and after the corner is shown for 
patches whose centroid speed is modified by a factor of oi~^ . These values of A7 are calculated from a root mean 
square average over q of the difference in action for all vortex patches with the same initial action and energy (and 
consequently separation from the wall) as the patch shown in Figure [T] The difference A7 decreases exponentially 
as a is increased (corresponding to e being decreased), which is consistent with adiabatic behaviour. 



are canonical variables for the Hamiltonian systents. 

SufRciently far before and after the corner, the strain experienced by the patch, and hence the action, are 
constant. We define /_ to be the action before, and /+ to be the action afterwards. The existence of an 
adiabatic invariant implies that the difference in action. A/ = /+ — /_, should be constant to all orders in e. 
Figure [Sfa) shows the evolution of / for the same motion as shown in Figure [T] During the interaction with 
the corner region / varies, but as t ^ ±cx) it tends to constant values, and the difference in these values is 
small (A/ ^ /), which is consistent with the existence of an adiabatic invariant. 

To confirm that A/ has the constant to all orders in e behaviour associated with the existence of an 
adiabatic invariant, we must analyse the variation of A/ with e. The rate at which E varies, and hence e, 
is determined by the speed at which the vortex patch moves through the corner. Therefore, we consider 
the motion of vortex patches with f and 9 as before (Equation [7]), but with x and y modified by a factor of 
a^^ . By varying a, we can vary the separation of the time scales e, which will vary as a~^. The resulting 
variation of A/ with a is shown in Figure IH^b). There is an exponential decrease as a is increased, which 
confirms that A/ is constant to all orders in e. 



The phenomenon of adiabatic invariance will not be limited to this particular boundary shape, but 
rather can be expected to occur for any boundary shape along which the motion of vortex patches sat- 
isfies Tstrain 3> Tj-ot • Adiabatic invariance could be broken by a boundary with variations in shape on a scale 
much less than the dimensions of the vortex patch, since Tstrain will be greatly reduced for such a boundary. 
However, fluctuations in the straining flow due to short scale boundary variations can be expected to decay 
exponentially with distance from the boundary on the same scale as the boundary variations (cf. solutions 
of Laplace's equation with sinusoidal boundary conditions). Thus the straining flow generated by such 
boundary variations will be exponentially small for any vortex patch whose separation from the boundary 
is large compared to its radius. 

We conclude that there are two possible scenarios for vortex patches in the vicinity of a boundary (in the 
absence of a background flow) : either the separation of the vortex patch from the boundary is comparable 
to its radius, in which case the vortex patch will undergo significant (order one) deformations to its shape 
(Figure [3]); or the separation of the vortex patch from the boundary is large compared to its radius, in 
which case the initial deformation to the shape of the vortex patch is small, and, due to the existence of 
an adiabatic invariant, will remain small. In particular, every time the patch returns to a straight section 
of boundary, its separation from the boundary and the periodically varying evolution of its shape must be 
very nearly the same. 
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V. CHAOTIC BEHAVIOUR OF THE ELLIPTIC MODEL 



The elliptic model for the motion of a vortex patch in the vicinity of a boundary is four-dimensional, and 
a generic feature of such four-dimensional dynamical systems is chaotic behaviour. In the examples that we 
have considered so far, any chaotic behaviour has been suppressed either by the motion being integrable, 
or by the presence of an adiabatic invariant. To demonstrate that chaotic behaviour is possible under the 
elliptic model, we briefly consider the motion of a vortex patch trapped in a corner by a background straining 
flow. 

If a pure straining flow, with stream function 5* = l3xy, is added to the quarter plane geometry of the 
previous section, then the straining flow can trap vortex patches in the corner provided the sign of /3 is 
appropriately chosen. Figure shows the streamlines for this background flow, and Figure IH^b) shows 
representative paths of vortex patches approximated as point vortices, which are trapped in the corner. This 
was previously noted by Llewellyn SmithSS. The additional excess energy due to the straining flow is 

Hb = -pLoAx.y, ~ |p(r, 9): E,,, (19) 

where 

E. = /3(j_°i) (20) 

is the strain rate of the background flow. 

The resulting motion of the elliptic vortex patches can be visualised via a Poincare section. Two such 
Poincare sections are shown in Figure [TUl for the case /3 = —0.028a;. Each Poincare section visualises the 
motion of several vortex patches, each with the same excess energy H, and depicts the points at which the 
angle 6* of a vortex patch passed forwards through 6 = 0. There is no chaotic behaviour visible for the lower 
value of H shown in Figure [rUT a). However, for the larger value of H shown in Figure [TUT b). two signs of 
chaotic behaviour are clearly visible: resonances and a 'sea of chaos'^^. Those patches whose motions form 
the the sea of chaos still orbit the central stagnation point, but the path of the patch centroid can differ 
significantly between consecutive orbits. Since the patch energy is constant, there must be a significant 
change of patch shape associated with this change of centroid position. 

The motion shown in Figure [TUT b) is actually unphysical due to the intersection of some of the elliptic 
vortex patches with the boundary of the domain; and we were unable to find any values of /3 and H 
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FIG. 10. Poincare sections in the Xc-yc plane of forward intersections with 6 = Q. The strain rate is /3 = —0.028a; 
and two different values of H are considered, (a) H/{uj'^R'^) — —0.44. The Poincare section closely matches the 
earlier point vortex paths, (b) H/{upR*') = —0.27. The Poincare section reveals chaotic behaviour. Resonances have 
opened up in some of the previously circular contours, and a 'sea of chaos' has formed around the edge of the region. 



that exhibited such clear chaotic behaviour without also being unphysical. Nevertheless, these results 
are informative since they emphasise the significance of the suppression of chaotic behaviour by adiabatic 
invariance that was observed in the preceding section. 



VI. INTERACTION WITH AIM ISLAND 



We now consider the deformation of a vortex patch in the vicinity of a circular island. This analysis has 
particular relevance to the propagation of Meddies in the Atlantic, which have been observed to interact 
with, and sometimes be broken up by, underwater seamounts^. 



A. No background flow 



We begin with the motion of a vortex patch in the vicinity of a circular island with no background flow. 
For the case of a point vortex, motion around an island can still be solved via the method of images. A 
point vortex of strength k located at x requires an image of strength —k located at 

X7=X— , (21) 

where a is the radius of the circular island^^. A side effect of this image is that it induces a circulation of 
—K around the island. Such circulation is a constant of the motion, and, since we are ultimately interested 
in vortex patches advected towards the island from infinity, we require it to be zero. The condition of zero 
circulation can be enforced via the introduction of a second image, of strength k, at the centre of the island. 
The resulting stream function is 

K / |x — X/| \ >^ , { 



V'(x) = — In ' , ^' = — In 1 - . (22) 



The image system for an elliptical vortex patch of uniform vorticity consists of a non-elliptical region of 
non-uniform vorticity. Fortunately, the elliptic model only requires the strain due to the leading-order point 
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FIG. 11. Motion of a vortex patch, initially circular with centroid (— 2.5-R, 0), around an island of radius a — R. The 
motion is visualised by snapshots of the patch location at time intervals of 40/a; for the elliptic model (dashed line) 
and the full contour dynamics solution (solid line). The path of the centroid under the elliptic model is also shown 
(grey line). 



vortex image, and thus the Hamiltonian (excess energy in the region + > a^) is easily calculated as 



where the strain rate, calculated from Equation (|22p . is given by 



1 - 



Zxy y —X 

y'^—x^ 2xy 



(24) 



The elliptic model is valid provided the centres of the vortex patch and its image are well-separated compared 
to the patch radius, i.e. Sc — a^/sc ^ R where Sc = {x^ + y^Y^'^- As before, this separation does not have 
to be particularly large before the elliptic model becomes a good approximation. 

Since the island geometry has rotational symmetry, there is an additional conserved quantity, which is 
the angular impulse 

uj j x-xd^x ^ bj {A{x\ + y\) + q-j^^, + . (25) 

J A 

The existence of an additional conserved quantity implies that the motion is integrable under the elliptic 
model. An example of the evolution of an initially circular vortex patch around an island is shown in 
Figure 1111 Also shown is the full contour dynamics solution, which is in good agreement with the elliptic 
model for t < 80/aj , but which starts to develop some non-elliptic deformations at later times. (See 
Appendix 1^ for a derivation of the contour dynamics method for a circular island.) 

We again use the maximum value of r for an initially circular vortex patch, rmax, as a measure of the 
deformation caused by the island. Figure [T^] shows this quantity evaluated for island radii of a = i?/2, R and 
2i? alongside the previous result for a straight boundary (a = oo) for comparison. For a given separation 
from the boundary, the deformation due to an island is less than that due to a straight boundary, and 
decreases as the radius of the island is decreased; this is due to a reduction in the magnitude of the straining 
flow. The strain rate due to an island of radius a has an asymptotic scaling of I s\ ^ oo, which 

decays faster than the l/y^ decay associated with a straight boundary, and which decreases in magnitude 
as the island radius is decreased. 
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FIG. 12. Maximum deformation rmax of an initially circular patch of vorticity in terms of the initial separation of its 
centroid from the boundary, Sc(0) — a, as given by the elliptic model. Larger islands lead to more deformation, and, 
in the limit a — >■ oo the result for deformation above a straight boundary is recovered. Solutions with Sc{Q) — a < R 
are unphysical since the vortex patch overlaps with the island. 



B. Background flow 

The effect of a background flow past the island can be investigated using an irrotational flow with the 
following stream function: 

n^.y)^Uy{l--^y (26) 

This represents a flow which tends to uniform flow away from the island. The associated streamlines 
are shown in Figure I13f a) , and a representative example of the resulting point vortex paths is shown in 
Figure [T3lb). Point vortices either pass above the island, pass below the island, or are trapped circulating 
around the island. However, such point vortex paths don't give us much information about the deformation 
of flnite area vortex patches, for which we need the elliptic model. 

The additional excess energy due to the background flow (using Equation and neglecting any moments 
higher than second order) is 



where 



Hk = -LoAUy, [ 1 - ^— ^ - -p(r, 9): Eb(xe, yj, (27) 

~^ Vc J ^ 

V) = (^2 + y2)3 (^Sx^y - y^ 3y^x - x^ ) 

is the strain rate at the patch centroid due to the background flow. 

Under the elliptic model, we can investigate the deformation of initially circular vortex patches released 
at a height j/g far upstream of the island [xg <^ —a). Depending on the strength of the background flow 
and the release height, the patches undergo one of three different motions: they either pass above the 
island, pass below the island, or pass sufficiently close to the island that the patch is indefinitely extended. 
Sufficiently close to the point vortex separatrix, another family of solutions may exist: vortex patches that 
loop around the island a number of times before continuing downstream. However, we were unable to find 
such solutions numerically, suggesting that any such solutions are limited to a very small range of release 
heights 2/0- Examples of the three different motions that were observed are shown in Figure [TH 
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FIG. 13. (a) Streamlines for the irrotational background flow '^{x,y) — Uy (l — a^{x^ + y^)~^) around an island 
of size a = R. (b) Paths of patch centroid under point vortex approximation with a = R and U/{u)R) = 0.1. The 
separatrix (heavy line) splits the motion into three regions: vortices passing below the island, vortices passing above 
the island, and vortices orbitting the island. This separation into three regions is generic and does not depend on 
the strength of the vortex. 

The deformation of an initially circular vortex patch can be investigated quantitatively by calculating the 
maximum deformation rmax as the patch passes the island. For a given ratio of island radius to patch radius 
(a/R), the elliptic model can then be used to quickly sweep over a range of ([/, j/o) parameter space. Such 
an analysis is shown in Figure [15] for the case a — R. Other island sizes give qualitatively similar results. 

The deformation of a patch depends significantly on whether it passes above or below the island. This is a 
consequence of the location of the point vortex separatrix (Figure [T3b), which passes closer to the top of the 
island than it does the bottom (for U /lj > 0). The path of a point vortex (of appropriate strength) gives a 
leading-order approximation to the path of the centre of the elliptic vortex patch. Thus patches starting just 
above the point vortex separatrix will pass closer to the island, and be deformed more, than those starting 
just below. The smaller the background flow speed, the further the point vortex separatrix lies from the 
island. Consequently, for small background flow speeds {\U\ ^ |w|i?), there is very little deformation of the 
vortex patch, independent of release height. Meanwhile, for large background flow speeds {\U\ \ui\R), the 
vorticity of the patch has very little effect on the motion and the centres of the vortex patches follow paths 
closer to those shown in Figure [T3la). In this 'ballistic' regime, those vortex patches that get deformed 
significantly correspond to vortex patches whose impact height is in, or close to, the range of impact heights 
covered by the island. 



C. Comparison with full solution 

The elliptic model is useful in allowing us to rapidly investigate large regions of (j/qi U) parameter space, 
but it is only an approximation to the full solution. For vortex patches that pass very close to the island, 
such as those that the elliptic model predicts are indefinitely extended, further investigation via calculation 
of the full solution is needed. 

The full solution, calculated for selected values of j/o with U/ (ioR) = 0.1 and a = i?, is shown in Figure [TBI 
From this figure, we note that the elliptic model is in good agreement with the full solution when r,jiax ^ 1-5- 
For those values of yo for which patches are indefinitely extended under the elliptic model (—2. Hi? < yo < 
0.61i? for this choice of parameters), the full solutions shown in Figure [TBI reveal two possible behaviours: 
patches are either distorted into a non-elliptical shape, before returning to an approximately circular shape 
downstream; or, for an intermediate range of release heights, they are broken up as they get extended around 
both the top and bottom of the island. 
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FIG. 14. Motion of initially circular vortex patches past an island of the same size (R = a) calculated using the 
elliptic model. The irrotational background flow has strength U/{luR) — 0.1. The motion of vortex patches is shown 
for three different values of the release height yo- Each motion is visualised by the path of the ellipse centroid 
alongside snapshots of the elliptic shape at time intervals of lO/cj. 



VII. CONCLUSION 

We have investigated the deformation of two-dimensional patches of uniform vorticity in the vicinity of 
fluid boundaries both via the Hamiltonian elhptic model of Melander, Zabusky and Styczek, and via contour 
dynamics simulations. For a straight boundary, the distance of the vortex patch centroid from the boundary 
is fixed and the motion is integrable under the elliptic model. In fact, the problem reduces to the well known 
problem of an elliptic vortex patch in a constant straining flow as analysed by Kidapi^. The elliptic model is 
a good approximation to the full solution even when the separation of the elliptic patch from the boundary 
is comparable to the dimensions of the patch. As the patch is moved away from the boundary, the straining 
flow generated by its image falls off as the inverse square of its distance from the boundary. Consequently, 
there is a rapid decrease in the deformation of the patch with distance from the boundary. 

The motion of a vortex patch in a quarter plane demonstrates nonintegrable behaviour under the elliptic 
model. However, due to a separation in time scales between the time scale on which the patch rotates 
and the time scale on which the straining flow experienced by the patch varies (the former being shorter 
than the latter), there is an adiabatic invariant. The existence of an adiabatic invariant acts to limit the 
deformation of vortex patches by preventing the accumulation of small deformations over time. We expect 
the motion of a vortex patch along boundaries with more complicated shapes to be similarly limited by 
the presence of an adiabatic invariant. When a background straining flow is added to the quarter plane 
geometry, vortex patches can become trapped in the corner region. Poincare sections of vortex patches in 
this system, evolving under the elliptic model, reveal the model's chaotic nature. 

We also investigated the deformation of a vortex patch in the vicinity of a circular island both with and 
without a background flow. Without a background flow, the motion of a vortex patch around a circular 
island is integrable under the elliptic model. The amount of deformation caused by an island is less than that 
caused by a straight boundary to a vortex patch at a similar separation, and increasingly so as the island 
radius is decreased. The amount of deformation also falls off more rapidly with distance from the island: the 
straining flow due to an island falls off asymptotically as the inverse fourth power of the separation between 
the patch and the island compared to an inverse square fall-off for a straight boundary. When a background 
flow is added, significant deformation of vortex patches is possible as they are driven towards the island 
boundary. The addition of an irrotational background flow that tends to uniform flow away from the island 
leads to three possible evolutions of vortex patches under the elliptic model: as they are advected towards 
the island, vortex patches either pass above the island, pass below the island, or are indeflnitely extended. 



17 



4 



2 - 




0.05 0.1 0.15 0.2 



U/{(oR) 

FIG. 15. Maximum deformation rmax of elliptic vortex patches with the same size as the island (R = a) over a range 
of background flow strengths U and release heights j/o- The region of parameter space for which vortex patches are 
indefinitely extended is shaded in dark grey. Elsewhere, contours of rmax are shown. For further information, the 
location of the separatrix for point vortices is shown (heavy line) as is the location of the island (light grey shaded 
region) . 



For those vortex patches that the elliptic model predicts are indefinitely extended, a closer examination using 
the full contour dynamics solution reveals two possible behaviours: either vortex patches are deformed into a 
strongly non-elliptical shape before passing around one side of the island and returning to an approximately 
circular shape downstream; or part of the vortex patch is split as part passes above the island and part 
passes below. 



Appendix A: Contour dynamics for a circular island 



Since the image vorticity distribution for a uniform vorticity patch outside a circular island has non- 
uniform vorticity, the derivation of the contour dynamics method for such a domain is more complicated. 
In fact, it is possible to apply the contour dynamics method to any finitely connected planar domain, as 
shown by Crowdy and Surana^, but here we just give a brief derivation of the result for a circular island. 

The streamfunction at x due to both the self vorticity and the image vorticity of a vortex patch of vorticity 
uj occupying an area A outside an island of radius a is given by 



|x-y||x| 



(Al) 



The velocity is given by derivatives of the stream function, and the contour dynamics method consists of 
converting the resulting area integrals into integrals about the vortex patch boundary. Taking a spatial 
derivative, and splitting into to contributions from the self and image vorticity gives 



Vx^(x) = -|-Vx^ln(|x-y|) d'y - Jjn ^ 



(A2) 



In order to rewrite the second term, which represents the image vorticity, as a boundary integral, we rearrange 
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FIG. 16. Motion of vortex patches with varying release heights driven towards an island by background flow; each 
panel shows a lower release height than the one before. The island and the patch have the same area, and the 
background flow has strength U/{ujR) = 0.1. Solutions were calculated using the elliptic model (dashed lines) and 
using contour dynamics (solid lines). The patch motion is visualised by snapshots of the patch position at intervals 
of 10/u;. In (a), the patch passes above the island and the elliptic model (rmax = 1.53) is in good agreement with the 
full solution; there is similar agreement in (f) when the patch passes below the island (rmax = 1.30). For intermediate 
release heights (b-e), the elliptic model predicts indefinite extension, as shown in (b). For these intermediate heights, 
the patch is either deformed into a non-elliptical shape while passing the island before returning to an approximately 
circular shape downstream (b,c,e); or the patch is extended around both sides of the island (d). 



it as 



Vx^(x) = -^Vx ^ In (|x - y|) <fy - Jjn^ 



|y| 



(A3) 



y - 

This allows both derivatives with respect to x to be rewritten as derivatives with respect to the dummy 
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variable y 



V.^(x) . V, m (|x - y|) d^y - |- V. (x^) • V, In 



xM 



d^y. 



(A4) 



Then, applying the divergence theorem, these area integrals can be rewritten as integrals around the bound- 
ary of the vortex patch 



Vx-^lx) = ^ / ln(|x-y|)fin 



UJ a 
2^W 



1-2- 



In 



OA 



y -^1 



dn. 



(A5) 



which is the desired contour dynamics formulation. 
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